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ABSTRACT 

We investigate the saturation level of hydromagnetic turbulen ce driven by the magnetorotational in stability 
in the case of vanishing net flux. Motivated by a recent paper of lBodo. Cattaneo. Mignone. & Rossi we here 
focus on the case of a non-isothermal equation of state with constant thermal diffusivity. The central aim of the 
paper is to complement the previous result with closure parameters for mean-field dynamo models, and to test 
the hypothesis that the dynamo is affected by the mode of heat transport. We perform computer simulations of 
local shearing-box models of stratified accretion disks with approximate treatment of radiative heat transport, 
which is modeled via thermal conduction. We study the effect of varying the (constant) thermal diffusivity, 
and apply different vertical boundary conditions. In the case of impenetrable vertical boundaries, we confirm 
the transition from mainly conductive to mainly convective vertical heat transport below a critical thermal 
diffusivity. This transition is however much less dramatic when more natural outflow boundary conditions are 
applied. Similarly, the enhancement of magnetic activity in this case is less pronounced. Nevertheless, heating 
via turbulent dissipation determines the thermodynamic structure of accretion disks, and clearly affects the 
properties of the related dynamo. This effect may however have been overestimated in previous work, and a 
careful study of the role played by boundaries will be required. 

Subject headings: accretion, accretion disks - dynamo - magnetohydrodynamics (MHD) - turbulence 



1. INTRODUCTION 

There are few concepts in classical physics that are 
equally fundamental as the conservation of angular momen- 
tum. One formidable consequence of this law is the for- 
mation of gaseous accretion disks around a wide range of 
astrophysical objects. Yet to explain observed luminosi- 
ties based on the release of gravitational binding energy 
dKing et alj 120071) . a robust mechanism is required to cir- 
cumvent the consequences of angular momentum conserva- 
tion. Being sufficiently ionised, these discs harbor dynami- 
cally important magnetic fields, which render the disk unsta- 
ble to a mechanism calle d magnetorotational instability (MRI, 
iBalbus & Hawlevll 19981) . releasing energy from the differen- 
tial rotation and converting it into turbulent motions. Ulti- 
mately, these motions are what is powering the redistribution 
of angular momentum on large scales. In the absence of ex- 
ternally imposed large-scale fields, a separate dynamo mecha- 
nism may be required to replenish s ufficiently cohe rent fields 
to drive sustained MRI turbulence (Vishniac 2009). This re- 
quirement becomes particularly apparent in models that ne- 
glect the vertical structure of the disk. In this case a con- 
vergence problem has been encountered ft ^essah et all 120071 
iFromang & Papaloizoul 120071: iBodo et aLl 120111) . which has 



been attributed to the l ack of an outer scale of the turbu- 
lence dDavis et al, 2010). An alternati v e expl anation has been 
suggested by [Kitchatinov & Riidigei] d2010h . who point out 
the problem of resolving the radial fine-structure related to 
non-axisymmetric MRI modes (required to circumvent Cowl- 
ing's "no dynamo" theorem). Convergence can naturally be 
recovere d when accounting for the vertical strati fication of 
the disk dShi et alj|2010t lOishi & Mac Lowll201 ll) . Together 
with rotational anisotropy, the introduced vertical inhomo- 
geneity can produc e a pseudo-scalar leading to a classical 
mean-field dynamo dBrandenburg et al1l!995al) . In a previous 
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paper (JGresselll2010l hereafter IG1 Oh . we have inferred mean- 
field closure para meters from isotherm al stratified MRI sim- 
ulations (also see Brandenburg 20 08|) and de monstrat ed that 
the "butterfly" pattern (see e.g. Miller & Stonel2000tlShi et alj 
120 lOt ISimon et all 120121) typical for stratified MRI can in 
fact be described in such a framework. In the current paper, 
we aim to extend this line of work towards a more realistic 
thermodynamic treatment, including heating from turbulent 
dissipation (Gardiner & Stone 2005; Piontek et al. 20Q9J) and 
crude heat transp ort. This new effort is largely initiated b y 
a recent paper of IBodo. Cattaneo. Mignone. & Rossil ( 12012b . 
hereafter iBCMRl The authors make the intriguing suggestion 
that the treatment of the disk thermodynamics will have a 
strong effect on the dynamo and accordingly on the satura- 
tion level of the turbulence . This is in contrast to the result of 
iBrandenburg et all (11995bl) . who find that turbulent transport 
is not affected by the presence of convection. 

As a first step to improve the realism of MRI simulations 
comp ared to the commonly applied isothermal equation of 
state, iBCMRl assume thermal conduction with constant ther- 
mal diffusivity, k. They then identify a critical value k c be- 
low which the primary mode of vertical heat transport changes 
from being predominantly conductive to being dominated by 
turbulent convective motions. One motivation of the present 
work is to scrutinize these models and at the same time obtain 
mean-field closure coefficients for the two suggested regimes, 
thereby confirming the assumption that the thermodynamic 
treatment of the disk affects underlying dynamo processes. 
A second focus of our paper will be the authors' choice of 
impenetrable vertical boundary conditions, which we find 
to have profound implications for the resulting vertical disk 
equilibrium. 

2. MODEL AND EQUATIONS 

The simulations presen ted in this paper extend the sim- 
ulation of iGress el (2010) to include a more realistic treat- 



ment o f thermodynamic processes as pioneered by Bo do et aTl 
d2012l) . Simulations are carried out using the second- 
order accurate Nirvana-iii code dZiegleri 120041) . which 
has been supplemented with the HLLD Riemann solver 
(Miyoshi & Kusano 2005) for improved accuracy. We solve 
the sta ndard MHD equations i n the shearing-box approxi- 
mation dGressel & Zieglerl2007l) employing the finite-volume 
im plementation of th e orbita l advection scheme as described 
in IStone & Gardiner! (120101); for interp olation we use the 
Fourier method by lJohansen et al.l(l2009l) . We here neglect ex- 
plicit viscous or resistive dissipation term s but include an arti- 
ficial mass diffusion term (as described in lGressel et al.l201lt) 
to circumvent time-step constraints due to low density regions 
in upper disk layers. We remark that, owing to the total- 
variation-diminishing (TVD) nature of our numerical scheme 
and the total energy formulation, heating via dissipation of 
kinetic and magnetic energy at small scales is accounted for 
even in the absence of explicit (or artificial) dissipation terms. 
Written in a Cartesian coordinate system (x, y, z) and with 
respect to conserved variables p, pv, and the total energy e = 
e + ^/ov 2 + |B 2 , with e being the thermal energy density, the 
equations read: 



d tP + v-(pv) = o, 



(i) 



dt(pv) + V-[pvv+p* - BB] =p [-V$ + aj] , (2) 

d t e + V- [(e + p*)v - (v-B)B] = pv- [-V$ + ai] 

+ V-(fcVT), (3) 



d t B 



Vx(vxB) = 0. 
V-B = 0. 



(4) 



with the total pressure given by p* = p + ^B 2 , and a fixed 
external potential $(z) = |il 2 z 2 . The inertial acceleration 
a; = 2S1 (qttx x-zxv) arises due to tidal and Coriolis forces 
in the local Hill system, rotating with a fixed $1 = ftoz, and 
where the shear-rate q = d In O/d In R has a value of —3/2 
for Keplerian rotation. 

We furthermore assume an adiabatic equation of state, such 
that the gas pressure relates to the thermal energy density as 
p = (7 — 1) e, with the ratio of specific heats 7 = I, as 
appropriate for a mono-atomic dilute gas. 

Finally, the temperature, T, appearing in the conductive en- 
ergy flux, is obtained via the ideal-gas law p = pT, where we 
chose units such that the factor /imH/ke relating t o the gas 
constant disappears. Following the approach taken bv lBCMRl 
we adopt a thermal conductivity, k, in terms of a constant dif- 
fusivity coefficient k, related via 



k = 



7 



7-1 



p K . 



(5) 



For the isothermal run, we do not evolve Equation ((3) but in- 
stead obtain the gas pressure via p = pT), with fixe d tempe ra- 
ture T = 1. Note that, in our units, T differs from lBCMRl by 
a factor of two, owing the alternative definition of the initial 
hydrostatic equilibrium, which is 



p(z) = p Q c 



2 /2H 2 



Poe 



-O 2 2 2 /T 



(6) 



in our case. For all simulations, we adopt the same box size 
of H x nH x 6H at a numerical resolution of 32 x 96 x 192 
grid cells in the radial, azimuthal, and vertical direction, re- 
spectively. This corresponds to a linear resolution of ~ 32 /H 
in all three space dimensions. 



TABLE 1 

Overview of simulation parameters, and results. 





EoS 


K [H 2 Q] 


vBC 


(M xy ) [10- a po] 


T m id [To] 


Ml 


isoth. 


- 


outfl. 


0.54 ±0.14 


- 


M2 


adiab. 


0.120 


outfl. 


0.63 ±0.17 


1.23 


Mi 


adiab. 


0.004 


outfl. 


0.81 ±0.21 


1.74 


M4 


adiab. 


0.004 


wall 


1.53 ±0.34 


4.94 



As typical for shearing box simulations, we initialize the 
velocity field with the equilibrium solution v = qil xy, and 
adiabatically perturb the density and pressure by a white-noise 
of 1% rms amplitude. The magnetic configuration is of the 
zero-net-flux (ZNF) type with a basic radial variation 



B = Bq sn\{2irx/ ' L x )z . 



(7) 



To obtain a uniform transition into turbulence, we further 
scale the vertical field with a factor (p(z)/p(0)) 1 ' 2 resulting 
in /3 P = const (= 1600 initially)]!] Owing to the divergence 
constrain, this of course can only be done by introducing a 
corresponding radial field at the same time. In practice, we 
specify a suitable vector potential. 

Horizontal boundary conditions (BCs) are of the stan- 
dard sheared-periodic type, and we correct the hydrody- 
namic fluxes to re tain the conservation prop erties of the finite- 
volume scheme dGressel & Zieglerl 120071) . For the verti- 
cal boundaries, we implement stress-free BCs (i.e., d z v x 



9 Z Vy 



0) with two different cases for the treatment of the 



vertical velocity component v z , namely: (i) impenetrable, and 
(ii) allowing for outflow (but preventing in-fall of material 
from outside the domain). We will demonstrate that this dis- 
tinction will have profound implications for the resulting den- 
sity and temperature profiles within the box. To counter-act 
the severe mass loss occurring in the case of open BCs, we 
continuously rescale the mass density, keeping the velocity 
and thermal energy density intactQ Unlike in earlier work 
dGressel et al.ll2012l) . which was adopting an isothermal equa- 
tion of state, we here do not restore towards the initial profile, 
but simply rescale the current profile. This is, of course, es- 
sential to allow for an evolution of the vertical disc structure, 
owing to heating via turbulent dissipation. We remark that 
such a replenishing of material can be thought of as a natural 
conseq uence of radial mass transport within a global disk. 

As in lBCMRL we use d z B = 0, B x = B y = as boundary 
condition for the magnetic field, and impose a constant tem- 
perature T = To at the top and bottom surfaces of the disk. 
The latter choice is motivated by the assumption that the upper 
disk layers are likely optically thin, and there exists a thermal 
equilibrium with their surroundings. As in previous work, we 
compute the density and thermal energy of the adjacent grid 
cells in the z direction to be in hydrostatic equilibrium. 

3. RESULTS 

The main motivation o f this pap er is to reproduce, as closely 
as possible, the results of lBCMRT where we then aim to estab- 
lish mean-field dynamo effects for the contrasting cases of ef- 
ficient versus inefficient thermal conduction, as studied there. 
Moreover, we shall begin to explore the impact of vertical 

1 Here and in the following, over-bars denote horizontal averaging. Addi- 
tional averaging in time is denoted by angle brackets, ( }. 

- Technically, this amounts to a volumetric cooling of the disk. 
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FIG. 1. — Time evolution of the volume-averaged Maxwell stress for the 
different models. For reference, we also list time averages in Table[T] 
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FIG. 2. — Time-averaged (for t 6 [50, 300] orbits) vertical profiles of the 
average Maxwell stress for the different models. 



boundary conditions, by studying the somewhat more realis- 
tic case allowing for a vertical outflow of material. 

We have performed in total four simulations: the two main 
simulations ('M2' and 'M3') adopt thermal diffusivity of k = 
0.12, and k = 0.004, respectively, representing the regimes 
of efficient and inefficient thermal transport (at molecular 
level), respectively. For both these simulations, and for a 
third isothermal reference run (model 'Ml' ), we adopt outflow 
boundary conditions. To assess the impact of the imposed ver- 
tical BCs (see column "vBC" in Table [T]), and to make direct 
contact with previous work, we adopt a fourth model, 'M4', 
with a value k = 0.004, and impenetrable boundaries at the 
top and bottom of the domain (labeled "wall" in the follow- 
ing). 

3.1. Comparison with \BCMIj 

We ran the different m odels for approximately 300 orbital 
times, 27rf2~ 1 ; note that IBCMRl use time units of O -1 in- 
stead. All time averages are taken in the interval t = [50, 300] . 
For the isothermal reference run, we obtain a time-averaged 
Maxwell stress, 



(Af a 



(B x —B x ) (By 



By)) 



(8) 



of (0.54±0.14) x 10~ 2 , which is comparable to the k = 0.12 
case with (0.63 ± 0.17) x 10~ 2 . In the case k = 0.04 
with outflow boundaries we find a somewhat higher value of 
(0.81 ± 0.21) x 10~ 2 , which is however significantly lower 
than in the otherwise identical case with impenetrable bound- 
aries with (1.53 ± 0.34) x 10~ 2 . While stresses are in- 
deed increased (by approximately 30%) at lower conductivity, 
clearly, the effect of the treatment of the boundaries is much 
more significant. 

Due to the strong fluctuations, the relative amplitudes are 
best seen in Figure [2] where we plot time-averaged vertical 
profiles of M xy normalized to the initial midplane gas pres- 
sure, Pq. Compared to the isothermal case Ml, the thermally 
conductive model M2 has a very similar vertical structure. 
The low conductivity model M3, where turbulent overturning 
motions dominate, is not too different in terms of its vertical 
profile either. Markedly, model M4, with impenetrable verti- 
cal BCs, shows maxima of AI xy around \z\ = 1.5 H . This is 
similar to the profiles shown in figure 8 of IBCMRl but note 
that their model shows a strong peak within 0.5 H of the ver- 
tical boundaries, while we only observe a very thin boundary 
layer in our model. 
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FIG. 3. — Space-time "butterfly" diagrams of the azimuthal magnetic field 
By for: model M4 with n = 0.004 and impenetrable boundaries (top), 
model M3 with n = 0.004 and outflow boundaries (middle), and model M2 
with k = 0.12 and outflow (lower panel). 



At the end of their result section, IBCMRl point-out that the 
spatio-temporal behavior of the dynamo is remarkably differ- 
ent between the k = 0.12, and k = 0.004 cases (see their 
figure 9). If we compare to run M4 with impenetrable bound- 
aries (uppermost panel of Figur e |3), we indeed find a simi- 
larly irregular butterfly diagram. IBCMRl argue that they find 
"no evidence for cyclic activity or pattern propagation" in the 
conductive regime. In contrast to this, looking at Figure [3] 
it appears that even model M4 at some level shows a prop- 
agating dynamo wave. More importantly, models M2, and 
M3 (which presumably are in the conductive, and convective 
regimes, respectively) show extremely similar dynamo pat- 
terns and cycle frequency (middle and lower panels in Fig- 
ure |3). This again suggests that the vertical boundaries have 
a profound effect on the mechanism of vertical heat transport 
and, as a consequence, on the dynamo. 
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4. — Time-averaged (for t g [50, 300] orbits) temperature profiles. 
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FIG. 5. — Same as Figure|4] but for the average density. 

To establish the fact that our model M4 is ind eed compara- 
ble to the corresponding simulation of lBCMRl we now look 
at ver tical pro files of the gas density and temperature. For the 
latter, BCMR had found a peculiar "tent" shape (see their fig- 
ure 4), i.e. a linear dependence of T on z, joined together at 
z = by a parabolic segment. In Figure |4] we compile the 
T(z) profiles, and such a tent-like profile can indeed be seen 
for model M4. However, in the case of low thermal diffusiv- 
ity, k = 0.004, we observe a much weaker deviation from the 
isothermal profile with outflow boundary conditions. 

Along with the tent-shaped temperature profile, iBCMRl 
found that heat transport from convection would erase the 
vertical density stratification and lead to a constant density 
near the disk midplane (see their figure 5). We observe a very 
similar density profile for model M4, which moreover shows 
strong density peaks at the domain boundary. These peaks 
are a consequence of enforcing the hydrodynamic fluxes to be 
zero at the domain boundaries. Unlike for impenetrable ver- 
tical boundaries, model M3 shows a more regular Gaussian 
density profile. For models M1-M3, the width of the bell- 
shaped density profiles is consistent with the trend in temper- 
ature and reflects hydrostatic equilibrium. Apparently, such 
an equilibrium cannot be obtained if solid-wall boundaries are 
a pplied. 

IBCMRl conjecture that when re crosses a critical value of 
k c ~ 0.02, the vertical heat transport changes from being 
mainly conductive to being predominantly convective. This 
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FIG. 6. — Conductive heat flux ( Fq ) (solid lines), and turbulent convective 
heat flux ( Ft ) (dashed lines) for the three non-isothermal models. 

was illustrated by their figure 6, where they compared the con- 
ductive heat flux 



F a = - 



7 _dT 

7 — 1 dz 



with the mean turbulent heat flux 



Fn 



7 



7-1 



pv,(T-T), 



(9) 



(10) 



where horizontal lines as usual indicate averages over the x 
and y directions. In Figure [6] we show the corresponding 
quantities for our non-isothermal runs: as in pr evious p lots, 
model M4 agrees very well with the result of IBCMRl but 
again differs significantly from the otherwise identical model 
M3 with outflow boundary conditions. Compared to M4, the 
net transport of heat is much reduced in the case of models 
M2, and M3. For clarity we plot these curves separately in the 
inset of Figure [6] with magnified ordinate. Much as expected, 
model M2 is dominated by a positive (i.e. outward) conduc- 
tive heat flux, whereas the turbulent flux is negligible. Unlike 
model M4 with solid-wall boundaries, the low-conductivity 
model, M3, only shows a very moderate level of convective 
heat flux, and notably one of the opposite sign. This negative 
flux appears to largely balance its positive conductive coun- 
terpart, implying low levels of net-conductive heat transport 
towards the disk surface and net-convective transport towards 
the midplane. It is instructive to note that the heating predom- 
inantly occurs at z ~ ±2 H, where the velocity and magnetic 
field fluctuations are highest, and not, as one might think, near 
the midplane. 

3.2. Mean-field coefficients 

One objective of this work was to test the dependence 
of the dynamo on the mechanism of vertical heat transport. 
One way of doing this is to establish mean-field closure pa- 
rameters for the new class of models with constant thermal 
diffusivity. L ike in previous work, we u tilize the test-field 
(TF) method (ISchrinner et al .1120051 120071) to measure coef- 
ficients such as the a effect, turbulent pumping, and eddy 
diffusivity. The used method (Brandenburg 2005) is "quasi- 
kinematic" dRheinhardt & Brandenburg 2010) in the sense 
that it has been found to remain valid into the non-kinematic 
regime in the absence o f magnetic background fluctuations 
dBrandenburg et al.l2008l) - whether this covers MRI is a topic 
of discussion. We here only briefly recapitulate the general 



framework of mean-field MHD and refer the reader to KjIOI 
for a more detailed descriptio n. A recent review about m ean- 
field dynamos can be found in lBrandenburg et ail (12012b . 

For the shearing-box approximation, due to its periodic 
character in the horizontal direction, there are no character- 
istic gradients expected in the radial or azimuthal directions. 
The natural mean-fields are accordingly those, which only 
vary in the vertical direction. With respect to the velocity 
u = v qflxy, the mean-field induction equation reads 



0.010 



d t B(z) = Vx [ u(z) x B(z) + £{z) 
+ (qQxy)xB(z)] , 



(11) 



where we have ignored a contribution due to microscopic 
magnetic diffusivity. Note that the explicit x dependence in 
the shear term, qilxy, drops out once the curl operation is 
applied. Furthermore, B z = B z (t = 0) = because of 
flux conservation in the periodic box. In this description, tur- 
bulence effects due to correlated velocity and magnetic field 
fluctuations are embodied in the mean electromotive force 



50) = u'xB', (12) 

which is typically parametrized as (lBrandenburgl2005l) : 

£i(z) = ocij(z) Bj(z) -fjij(z) £jki d k Bi(z) , 

where?, j G {x,y}, k = z. (13) 

Given explicit knowledge of the rank-two a and fj tensors, this 
closure allows to formulate the mean-field induction equation 
(TTTT) in terms of mean quantities alone, leading to the classical 
a£l dynamo descr i ption, fi rst applied to MRI turbulence by 
iBrandenburg et all (fl995al) . 

Tensor coefficients representing the closure parameters are 
presented in Figure [TJ for model M2, with a high value k = 
0.12 of the thermal diffusivity. The obtained profiles are 
largely similar to the corresponding curves of the isothermal 
model Ml (not shown), which more over agreqj with previous 
isothermal results reported in iGlOi Notably, for model M2, 
the contrast between the disc midplane and the upper disc 
layers is somewhat less pronounced compared to the purely 
isothermal case. This trend is continued when going to lower 
therm al diffusivity (see Figure [8] below). Unlike reported in 
IGIOL we now find a xx and a yy to be predominantly of the 
same sign, which would argue in favor of a kinematic (rather 
than magnetic) origin of the effect. Note however the signif- 
icant fluctuations in a xx , which cast some doubt on whether 
this coefficient can be meaningfully determined in the pres- 
ence of shear. On the other hand, a fluctuating a should not 
be disregarded as a possible so urce of a mean-field dynamo 
(IVishniac & Brandenburgl 19971) . 

Before we proceed, we briefly discuss the remaining coef- 
ficients. In panel (b) of Figure [7] we show the off-diagonal 
tensor elements of the a tensor, which are dominantly sym- 
metric, i.e. a xy ~ a yx . We remark that for the clas- 
sical diamagnetic pumping effect, one would require anti- 
symmetric parity. The observed symmetry may however be 
interpreted as differential pumping, i.e. transporting radial 
and azimuthal field in opposite directions. For reference, we 
plot the mean vertical velocity u z (see dashed line), which 
additionally transports the mean field and hence leads to the 

3 This is with the exception of a xx , which generally appears to be poorly 
constrained by our simulations. 
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FIG. 7. — Mean-field coefficients computed via the TF method for model 
M2 with k = 0.120. Axis labels indicate curves plotted in dark (a xx , . . . ) 
or light (ctyy, . . . ) colours, respectively. The mean vertical velocity, u z , in 
panel (b), and the classical estimate for the turbulent diffusion, in panel (c), 
are shown as dashed lines. 



characteristic acceleration in the butterfly diagram. The diag- 
onal parts of the fj tensor are shown in panel (c), where we also 
plot the rms velocity fluctuation (dashed line). Apart from the 
boundary layers, the turbulent diffusivity agrees well with the 
theoretical expectation 



77 T 



(14) 



that is, assuming a coherence time r = 0.03 il^ 1 of the turbu- 
lence. Given the dominance of the azimuthal field, the coeffi- 
ci ents fj xx , and f) vv are surprisingly iso tropic (unlike predicted 
bv lVainshte in. Parke rT& Rosnerll993l) . It is however interest- 
ing to note that while fj xy , shown in panel (d) of Figure UJ is 
identical to the diagonal elements of the diffusivity tensor in 
panel (c), its counterpart ij yx is much smaller. With negative 
shear and both coefficients positive, the dynamo based on the 
Radler (1969) e ffect i s decaying (cf. the dispersion relation 
in Brandenburg 2005, appendix B). This does however not 
exclude the possibility that fj yx has an effect on the overall 
pattern propagation. 

We now proceed to the corresponding coefficients for the 
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FIG. 8. — Same as Figure [jj but for the model M3 with k = 0.004. Note 
that for the sake of direct comparison, axis ranges are kept fixed with respect 
to the previous figure. 




FIG. 9. — Same as Figure [7j but for the model M4 with k = 0.004 and 
impenetrable vertical boundaries. Note the different axis ranges as compared 
to the previous two figures. 



case of inefficient thermal conduction. In Figure[8]we accord- 
ingly show the a and fj tensor components for model M3. 
We recall that in the quasi-isothermal case (cf. Figure |7) the 
a xx and a yy coefficients showed a trend to flatten and even 
reverse their slope near the midplane (also cf. figure 9 in 
lBrandenburgl2008l) . This was reasoned to be related to a neg- 
ative a effect due to m agnetic buoyancy dBrandenburgll 1 99a : 
iRiidiger & Pipinll2000h . In contrast, here the a xx curve shows 
a more monotonic dependence on z, indicating that such mag- 
netic effects may be less pronounced in this case. Such a 
trend appears consistent with reduced magnetic buoyancy in 
the case of a stiffer effective equation of state. Moving on to 
panel (b), we note that a xy is now suppressed and even shows 
a slight trend to change its sign - indicating a possible signif- 
icance of diamagnetic pumping in the adiabatic case. For the 
turbulent diffusivity plotted in panel (c), we observe a signifi- 
cant deviation from the classical expectation (dashed line), re- 
sulting in a nearly constant t]t{z)- We conclude that the equa- 
tion of state and the means by which energy is transported to 
the upper disk layers indeed have subtle effects on the inferred 
dynamo tensors. Which of the differences seen between Fig- 
ures [7] and [8] is in the end responsible for the enhanced dy- 



namo activity seen in model M3, will require further careful 
study. For completeness, in Figure [9] we also show the re- 
sults for model M4, where dynamo coefficients are larger by 
a factor of several. This can be considered consistent with the 
much higher Maxwell stresses observed in this case. Unlike 
for model M3, fj(z) now shows a pronounced z dependence, 
and a yy , an d a vx are in fact ne gativJl (for \z\ < 2 H) as 
suggested by iBrandenburd (1 1 998b . 

4. CONCLUSIONS 

The primary goal of o ur work was to make contact with re- 
cent results by [BCMR J and to complement their work with 
a direct measurement of mean-field dynamo effects via the 
TF method. Given that we have used a very similar numer- 
ical method and applied identical parameters (i.e. for model 
M4), it should not surprise the reader that we can satisfac- 
torily co nfirm all aspects of the corresponding simulation by 
iBCMRl . Minor discrepancies arise with respect to the bound- 
ary layers, which may be related to the detailed treatment of 
the hydrostatic equilibrium there. 

4 i.e. for z > 0, and accordingly positive for z < 



Accordingly, we can confirm their main result, namely that 
- in the presence of impenetrable vertical boundary condi- 
tions - one observes a transition from a conductively dom- 
inated vertical heat transport to a state that is regulated by 
convective overturning motions. This transition obviously de- 
pends on the value o f the app lied constant thermal diffusivity 
k. Like reported in BCMR, in the k = 0.004 case, we ob- 
serve a flat density profile (even with a slight minimum at the 
disk midplane), and a "tent"-shaped temperature profile - pre- 
sumably established by convective heat transport as a result of 
Rayleigh-Taylor-type instability. This case is also associated 
with a much increased dynamo activity (by an order of mag- 
nitude in the TF coefficients, not shown here), resulting in an 
overall Maxwell stress that is increased by a factor of three 
compared t o the isothermal reference model. Unlike specu- 
lated by BCMR, we however do not think that the enhanced 
dynamo activity is related to the magnetic boundary condi- 
tions. This is despite the fact that such a connection indeed 
exists for the unstratified case (K apvla & Korpill201 ll) . where 
the different magnetic boundary conditions serve to create an 
inhomogeneity in an otherwise translationally symmetric sys- 
tem. We rather attribute the different dynamo regime to the 
overall different hydrodynamic state - which however appears 
largely influenced by the choice of impenetrable boundary 
conditions. 

A separate set of models (M1-M3) with a more natural con- 
dition allowing the gas to flow out of the domain shows much 
less dramatic effects when going to the low thermal diffusiv- 
ity regime. Naturally, one arrives at moderately hotter disk 
interiors along with more spread-out, yet still Gaussian den- 
sity profiles. Dynamo TF coefficients are somewhat altered in 
this case, along with a roug hly 30% higher tu rbulent Maxwell 
stress. Establishing a link (lBlackmanll20 1 Oh between the tur- 
bulent transport coefficients in the momentum equation (i.e. 
the Maxwell and Reynolds stresses) and the induction equa- 
tion (i.e. the a effect, turbulent diffusion, etc.) will be key 
to understand magnetized accretion in a quantitative manner, 
and derive powerful closure models in the spirit of lOgilviel 



( 12003b or [Pessah et al. (2006). A possible direct extension 
of the existing models with varying amounts of thermal con- 
ductivity may be to cross-correlate (M xy ) with e.g. a yy for 
va rious values o f k. S uch a connection has been suggested 
by iBrandenburgl (119981) and been d erived in the quasi-linear 
regime by [Rudiger & Pipinl (|2000). A complication in this 
endeavor however arises from the fact that the TF coefficients 
are likely measured in a magnetically affected, i.e. quenched 
state. 

Given the dramatic effect of open versus closed vertical 
boundaries demonstrated in this paper, it will be of prime 
interest to study the connection betwee n the disk and the 
launching of a magnetically driven wind (lOgilviell2012b . in- 
cluding a possible influence of the wind on the d i sk dynamo . 
The amount of recent work ( From ang et al.ll20T2"l lMollll20 1 2l 
Bai & Stone 2012; Lesuret al .1 120 131) illustrates the impor- 
tance of this issue. In terms of the vertical boundaries imposed 
on the temperature, including a transition into an optically- 
thin disk corona will be important. Then a radiative boundary 
condition consistent with black-body radiation can be applied. 
To conclude, we want to emphasize that, clearly, the pre- 
sented simulations can only be regarded as a first step towards 
a realistic treatment of the disk thermodynamics. Ideally, 
full-blown radiative transfer should be employe d, and sim- 
ulatio ns of radiation dominated accretion disks (Bl aes et ail 
12011b demonstrate that this has indeed become feasible. 
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